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, Abstract. A model for the production of quarkonium states in the midrapidity 

'^J ■ region at RHIC and LHC energy range is presented which explores well understood 

properties of QCD only. An increase of the quarkonium hadronisation time with 
the initial energy leads to a gradual change of the most important phenomena from 
^ ' fixed target- to collider-energies. We evaluate nuclear effects in the quarkonium 

OO , production due to medium modification of the momentum distribution of the 

C 2 ^ . heavy quarks produced in the hard interactions, i.e. due to the broadening of the 

f^^ ' transverse momentum distribution. Other nuclear effects, i.e. nuclear shadowing 

^\ ' and parton energy loss, are also evaluated. 

o 
o, 

1. Introduction 

At fixed target energies the liadronisation time tf for QQ pair is small as compared 
y ■ to the nuclear radii, at collider energies this time is much larger than the size of the 

pH I nucleus. Thus, at fixed target energies, the interaction of heavy quarkonium with the 

surrounding medium can be described effectively as hadronic W p| even within the 
hadronisation time. However, at collider energies a description of the process in terms 
of hadron degrees of freedom becomes questionable since a bound state of the heavy 



X 



JH I quark pair is formed far outside of the nucleus. 

One major effect at collider energies is the transverse momentum broadening of 
QQ pairs travelling through a medium. This changes the overlap integral between 
the QQ pair and the wave functions of the heavy quarkonia. It is important that soft 
interactions of heavy quarks with nuclcons do not change the longitudinal momentum 
of heavy quarks which may appear sensitive to the interaction with comovers which 
are neglected in this paper. Nuclear effects in the momentum distribution within the 
heavy quark pair due to the radiation of gluons which is sensitive to the properties of 
the QCD-matter produced in a nucleus-nucleus collision at collider energies ||^ will be 
estimated separately in the section six. 

The typical time scales for the creation of a hidden heavy flavour state are 
investigated in section two, in the third section a model is developed to calculate 
the influence of nuclear effects on the production of states with hidden heavy flavor 
in the midrapidity region at collider energies. In the fourth section this model will be 
used to predict heavy quarkonium products in AB collisions at RHIC and LHC and 
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in section five and six quarkonium suppression due to nuclear effects on the parton 
distribution functions and parton energy loss are evaluated. 

2. Coherence Lengths 

There are different time scales relevant for the production of quarkonium states: 

1) the time needed to produce a heavy quark pair in a hard collision, 

2) the time needed for a QQ pair to form a bound state. 

The production time of a QQ pair in its rest frame is given by Tp = ^— . This is 

0.13 fm/c for cc and 0.05 fm/c for bb. The Lorentz factor of the pair at midrapidity 
in the rest frame of the target is 7 = -^^ w 10, 20, 100 and 3000 for a Quarkonium 
state in the midrapidity region at SPS, Fermilab, RHIC and LHC energies. Thus, at 
SPS and Fermilab fixed target energies 'ycTp is smaller than the average internucleon 
distance in nuclei r^N ~ 1-8 fm. Thus, the production of heavy quark pairs is 
incoherent. However, at RHIC the production distance of cc pairs is already as large 
as the diameter of a gold nucleus, and for bb pairs CTp > 1.8 fm, but this is still small 
as compared to the nuclear radius. At LHC both production distances exceed the 
diameter of a lead nucleus by an order of magnitude. 

The hadronisation time t^ resp. the coherence length l^ of heavy Quarkonium is 
ic = c- t H = ^^^ with: 

AE ^ ^p^ + (M,, + AMy - ^7^%-^ ^ ^-00^-f-Ka ^ M^ ^ ML . 
Here p is the momentum of the Quarkonium in the rest frame of the target, Mqq is 
the mass of the QQ-pair and AM = J ip'^{k)^d^k/ J ilj'^{k)d^k , where tp{k) is the 
wavefunction of the Quarkonium state in momentum space. That means AM is the 
average kinetic energy of the QQ-pair in the bound state and Ic/'y = 0.44(0.34) fm for 
the J/* (T). 

Thus, for charm and bottom production at RHIC and LHC Ic > 2- Ra {Ra is the 
nuclear radius), Zc < 2 • Ra for fixed target at Fermilab (7 ~ 20) and at SPS energies. 
The applicability of the approach developed in this paper requires that ^c > 2 • Ra 
which is fulfilled at RHIC and LHC (and at large Xf at Fermilab fixed target energy 
range, where 7 ~ 100 also can be reached). 

3. Description of the Model Approach 

We assume here that QQ pairs are produced in AB collisions predominantly in hard 
collisions. The basic quantity is the cross section of production of QQ pair with 
light cone momenta Zi,ki, which we parametrize as d2fc dz^dz d^k ~ Dab{zi, Z2) ■ 
ex-p{—B{AB){kl + A:|)). Here ki {i — 1, 2} are the transverse momenta of the Q and 
the Q quark and z^ {i = 1, 2} are the fractions of their longitudinal momenta. Such a 
factorization does not contradict the data in pp collisions M. 

To evaluate the suppression of hidden heavy fiavour production resulting from 
the broadening of the transverse momentum distributions of Q quarks due to final 
state interaction, we deduce first a relationship between the slopes for the various 
processes of heavy quark production. In the following we use the relative transverse 
momentum kt = ^^ ^ and the total transverse momentum pt = ki + k-z of the pair, 

writing 'Xf-rpZ^,fJ = Dab{zi,Z2) exp (^B{AB) (-f - 2k? 
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To take into account possible nuclear effects on the longitudinal momentum 
distribution we make the ansatz Dab{zi, Z2) = D{AB) ■ fAsiPz) -exp I — „/ ) . where 



Pz and kz are the total and relative longitudinal momentum. We further assume that 
fAsiPz) — fppiPz) , which means that we neglect parton energy losses of the pair, 
this effect will be discussed later on. The normalization condition follows from the 
QCD factorization theorem Q for the total cross section: bIaby^-c — ~ B(pp)'^^c 
The differential cross sections are proportional to the square of the two body wave 
functions cj). 

The production cross section of bound states of heavy quarks is proportional to 
the overlap integral of the two-body wave function and the wave function of the bound 
il){kt) state to get ^ 'y{AB^Quarkontum+x) ^ \{tp{kt,kz)\(j)AB{kt,Pt,kz,Pz))\^ . Here wc 
neglected the difference between the current quark mass in the two body wave function 
and the constituent quark mass in the wave function of the bound state. With this 
one can evaluate the survival probability: S = -^ — -, — ,"*"'"' ""'""' — _L^. Our final 

^ *^ AB-a{p+p — >Quarkonium+X) 

result is then 



^ ^ B{AB)Cab 



B{pp)C, 



pp 



J d^kiPik) exp(-B(AB)fc2) exp (- 



2C 



J d^k^{k) exp{-'B{pp)kf) exp ( - 



k^ 

2ch: 



(1) 



up to nuclear effects in the structure function. 

Note that if one defines the survival probability as the ratio of the differential cross 
section ^^ for nuclear and nucleon targets their pt dependence would be a factor 

exp ( — '-^ — ^^P( J. That means that the pt dependence of J/^ suppression is 

due to the broadening of the transverse momentum distribution as a result of the final 
sstate interactions of the Q quarks in the nuclear medium. 

4. Predictions for Collider Energies 

To adjust the parameters of eq. (II]) we use the quarkonium wave functions of ref. p, p| 
and the following data: 

The E789-group (Fermilab) measured a value of B{pAu) = 0.91 ± 0.12 GeV"^ 
with the fit function j^{D,D) ex exp{—B{pAu) ■ k^) in proton-gold collisions for 

the transverse momentum distribution of D and D mesons. The same group found 
n = 0.55 ± 0.01 GeV"^ @ with |^(J/*) ex exp(-n • p'^), also in proton-gold for the 
transverse momentum distribution of J/^I^'s. Within the error bars these values of n 
and B{pA) satisfy the relation n ~ ^^ ' that was predicted in the previous section. 
With the fit §1 to the pt broadening of J/*'s this yields B{pp) = 1.18(1.5)GeV"^. 

There are no data for the nuclear broadening of the longitudinal momentum 
distribution. We assumed isotropic production in the cc rest frame which yields 
C^p = 1/(2 • B{pp)) « 0.42 GeV^. In fig. 1 the result of eq. (|) is plotted versus 
the transverse momentum broadening of the J/'^'s: Ap"^ = 2/B{pp) — 2/B{AB). 
Ap( — 0.48 GeV^ was found at Fermilab energies in pAu. C\g = Cpp is used as a 
first approximation. S ~ I for the J/^P. That means there is practically no change 
due to the broadening of the transverse momentum distribution. 

For the T meson we use < p? >= 2.59 GeV^ |l§ and thus B{pp) = 0.77GeV~^ 
and Cp = 1/(2 • B{pp)) k, 0.65 GeV^. For the T meson production we obtain even a 
slight enhancement. 
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Figure 1: The righthandside of 
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A(Pt^)=PtAB^-Ptpp^ in GeV^ 

Ref. [gj reported that Ap^ is nearly a factor 5 larger for J/^'s than for Drell 

Yan pair production. In our model this is due to a color factor of | because cc 

production is dominated by gluon-gluon fusion, while Drell Yan is due to quark- 

antiquark annihilation and due to the factor of two between the transverse momentum 

distribution of c resp. c Quarks and the bound states evaluated above. 

5. Nuclear Effects on the Gluon Distribution Functions 

The term nuclear effects on the gluon distribution functions means that the ratio of 
the gluon distribution functions in a nucleus A and a nucleon TV differs from 1, i.e 
AGM^ifQ^) ^ ^ ^^ ^^^ shadowing region, ^G^lfgl^ > 1 for 0.1 < a; < 0.3 and for 
X > 0.8. The Bjorken x's of interest for the production of J/^'s x — 0.0155 and T's 
X = 0.0473 are x(QQ) = ^. 

Up to now, the transition points between these regions are not well known. 
Especially for gluons, the error bars of the data are large. Also from a theoretical 
point of view they are not unambigious. Therefore, the calculation of these nuclear 
effects, at not too small x, is model dependent. To reduce this model dependence we 
use here two different recent approaches from ref. O, O (called EKS in the following) 
and from ref. Il^ (FS). (Nuclear shadowing in (FS) has been calculated in a nearly 
model independent way, but the region of enhancement requires modelling.) 

In tab. 1 one can see that quarkonium production at RHIC is enhanced while 
it is shadowed at LHC. At LHC energies the predictions of both approaches become 
qualitative different. The latter one based on significantly larger cross section of 
diffraction in the gluon channel or equivalently on the 9/4 enhancement of color dipole 
build of color octets. It predicts a stronger shadowing of gluons than the shadowing 
of quarks and of antiquarks, while in the EKS model it is as strong as the shadowing 
for sea-quarks. 

6. Suppression due to Parton Energy Loss 

Parton energy loss shifts the partons to smaller x. We use here the parton energy loss 
model reviewed in ref. M. AE is the parton energy loss in a medium. For cold QCD- 

matter AE = -4 GcV • f — ^1 . For hot QCD-matter AE = -60 GeV • f — ^ 
Vio tmy ^ Vio im^ 

With L = Rau = 7 fm this yields AE = —2 GeV for cold matter and 
AE = —30 GeV for hot matter. For RHIC energies this value has to be divided 
by 100 since in Q the parton energy loss was calculated in the rest frame of the target 
while we are interested in the midrapidity region. Therefore, an initial parton with 
x{J/'^) = 0.017 in hot QCD-matter is needed to produce a J/^. For pA collisions 
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y/s in GeV 


60 


140 


200 


6000 


J/* (FS) 


1.11 


1.05 


1 


0.65 


J/* (EKS) 


1.16 


1.01 


0.96 


0.84 


T (FS) 


1.09 


1.04 


1.09 


0.83 


T (EKS) 


1.104 


1.138 


1.12 


0.87 



Table 1: 



initial energies (see text) 



at different 



this yields an suppression factor of -^^ — — — T-^-^rrz ~ 0.87 for the gluon distribution 

^ ^^ G(0. 0155,3.1 GeV) ^ 

parametrization of Ref. |14, n5|. For the T this factor is 0.94. For cold QCD-matter 
this is negligible, since x — 0.0156 for the J/vf. This yields a suppression factor of 0.99. 
For AB collisions this values have to be squared, since both initial gluons can loose 
energy. Therefore, the suppression for the J/^ (T) in hot matter in AB collisions is 
0.76 (0.88). This educated guess must suffice for the time being, because the parton 
energy loss is theoretical not well understood. 



7. Discussion and Conclusion 

• the relevant physics changes due to coherence length: 

— > hadronic description at fixed target, partonic at colliders 

• nuclear effects on the parton distribution functions yields: 
— > anti-shadowing at RHIC, shadowing at LHC 

• parton energy loss seems to be unefFective in cold QCD matter 

• other possible effects not discussed here include: 
— > interaction with comovers 

— > long living phase of hot QCD-matter 

— > J/^I'-production and T comes also from the decay of higher resonances 

• pA collisions are needed to distinguish between the different effects. 
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